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Abstract 

^ ' In this paper we develop an analytical model in order to study electromagnetic 

■ processes involving loosely bound neutron-rich and proton-rich nuclei. We con- 
Q i struct a model wave function, to describe loosely bound few-body systems, having 
0^ ' the correct behaviour both at large and small distances. The continuum states are 

■ approximated by regular Coulomb functions. As a test case we consider the two- 
, body Coulomb dissociation of and, the inverse, radiative capture reaction. The 

difference between using a pure two-body model and the results obtained when 
incorporating many-body effects, is investigated. We conclude that the interpreta- 
tion of experimental data is highly model dependent and stress the importance of 
^ ' measuring few-body channels. 

c ■ 



Key words: 

X ■ PAGS: 21.60.Gx; 25.40.Lw; 25.70.De; 27.20.+n 



Electromagnetic processes such as radiative capture, photo-dissociation and 
Coulomb dissociation have always been, and still remain, an excellent tool to 
investigate nuclear structure. Unfortunately, for radioactive nuclei the radia- 
tive capture experiments are very difficult and studies of photo-dissociation 
are virtually impossible. However, with the advent of radioactive beam facili- 
ties, nuclear structure of dripline nuclei can be studied using Coulomb dissoci- 
ation on heavy targets. The Coulomb breakup of these nuclei is of interest also 
in nuclear astrophysics, since it can be related to the corresponding radiative 
capture process at astrophysical energies [1]. In this paper, we will present 
an analytical model for Coulomb dissociation (and consequently for radiative 
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capture) of loosely boTind dripline miclei (see also our recent conference con- 
tribution [2]). We want to stress that our approach will be general in the sense 
that both neutron-rich, and proton-rich systems can be studied. The possi- 
bility to study reactions analytically is very appealing since it often allows a 
deeper physical understanding of the process. For some cases such studies are 
not only possible but might even give results that are directly comparable to 
experimental data and more advanced, numerical calculations. In particular, 
this is the case for electromagnetic transitions between a loosely bound nu- 
clear state and a pure Coulomb continuum. "Loosely bound" implies that the 
nucleus will exhibit a large degree of clusterization. We will consider the case 
of two clusters, but will also discuss the three-cluster channel (see also similar 
approaches for one-neutron [3,4,5] and two-neutron [6,7,8] halo nuclei). After 
a general description of our model we will exemplify it with an application to 
^B. Our main focus will be the nuclear structure effects, and we will utilize 
the advanced three-body model of Grigorenko et al. [9], while the Coulomb 
dissociation is considered only to first order. As it was shown recently [10,11], 
higher order effects should not be significant at beam energies higher than 
around 70 MeV / A. 

Our starting point for calculating electromagnetic cross sections will be the 
EA strength function for a transition from an initial state to a final, continuum 
state with energy E 

^^^-^J^T.J drf\{f\\M{EX)\\t)f6{Ef-E), (1) 

where dr/ is the phase space element for final states, M.(EX) is the electric 
multipole operator and \ f) are the initial and final states in the center of 
mass subsystem. 

We will consider loosely bound systems of two clusters {c + x) and, in particu- 
lar, we will study transitions to the low-energy continuum in which excitations 
are manifested as relative motion between the clusters E — fi^k'^ /2iicx, where 
IXcx is the reduced mass of the system. Introducing the intercluster distance, r, 
the corresponding cluster EA operator (operating only on the relative motion 
of clusters) is 

A1(EA,;u) =eZ(A)r^rAM(r), (2) 
with the effective multipole charge Z{\) = ij,^^{Zx/m^ + {—lyZc/m^). 

The strength function is the key to study several reactions. The photo-dissociation 
(^4 -|- 7 — >• c -|- a;) cross section is given by 

Kx.p. (27r)^(A + l) fE,^''- UB{EX) 



A[(2A + 1)!!]2 V^cy dE 
where the photon energy, E^ — E + Eq, is larger than the binding energy, Eq. 



2 



The inverse radiative capture reaction can be studied using detailed balance 



RA/m {^^\ 2{2Ja + 1) _EA^F^ (AX 

where J, is the spin of particle i. Using first order perturbation theory, and 
the method of virtual quanta [12,13], the energy spectrum for Coulomb dis- 
sociation on a high-Z target can be written as a sum over multipole, ttA, 
photo-dissociation cross sections multiplied by the corresponding spectra of 
virtual photons, n.j^x{E^), 

dac ^ nEA(^7) ^EA^F ^ I V^ ^MA(-g7) ^MA^F ^ 
EA -^7 MA -^7 

Note that since MA transitions are usually [13] strongly suppressed we will not 
consider them in this work. In order to study all these reactions analytically 
we will propose a model function to describe the (loosely) bound state of a 
two-body system. We will only be interested in direct transitions to a "clean" 
continuum, i.e., with all nuclear phase shifts equal to zero. Thus, the final 
state, with angular momentum 1^ between the clusters, will be described by a 
regular Coulomb function 

<i>if r) = \/f ^^''e'^'^ (k, r) , where (6) 
Fi^{k, r) = Ci^{k)e''''{kryf+\F,{lf + 1 + t7]{k); 2lf + 2; -2ikr), (7) 

and (Tij is the Coulomb phase, r]{k) — ZcZ^^ HcxjY^k is the Sommerfeld pa- 
rameter, Ci^{k) = 2'/e-'^''('^)/2|r(;^ ^ 1 ^ iri{k))\l{2l^ + 1)!, and ^F^{a\h\z) is 
the confluent hypergeometric function [14]. 

The reduced matrix element introduced in the definition of the strength func- 
tion, Eq. (1), contains a radial integral. With our approximation for the con- 
tinuum state this integral takes the form 

("OO 

Ii^{k) = J dre-'''W^f+\Fi{lf + l-ir){k)]2lf + 2]2ikr)r^(j)b{r). (8) 

Here, 06 (r) is the two-body, relative motion wave function (WF) describ- 
ing the initial state. For large r, and with angular momentum Zj between 
the clusters, this radial function should be proportional to the Whittaker 
function ;._|_i/2(2fi:or) (see Ref. [14]), where rjQ = ZcZ^e^Hcx/h^KQ and 

£^0 = 1^0 / 2 Hex is the binding energy. 

In most studies on loosely bound systems, the Whittaker function has been 
used to describe the bound state for all r. However, this approximation is only 
motivated if the transition matrix element is dominated by contributions from 
very large r. This should be the case for reactions at very small energies. For 
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real experimental energies {E > 100 keV), the WF of the bound state should 
be constructed in a more realistic way. Therefore we will introduce a model 
function that describes the bound state c + x WF accurately for all distances. 
This is done by considering the behavior at small and large r. We have already 
pointed out the the WF should be described by a Whittaker function at large 
r. Furthermore, for a two-body system, consisting of point-like particles, it 
should behave as r''+^ as r ^ 0. Both asymptotics are fulfilled using the 
following model function 



^:rf "'V) = ^W -,„,,+v2(2«:or) (l - e"'^-)''''' , (9) 



where A'^ is the normalization constant and 7 denotes the parameters {k.q, Ki, rjo}. 
The parameters Kq and 770 ^-re defined by the binding energy, charges and 
masses, while Ki can be fitted to give the correct size of the system. Using 
this WF, and solving the integral (8) numerically, it is possible to get very 
good estimates for the electromagnetic reaction cross sections. As a remark 
we want to point out that, for a one-neutron halo nucleus (770 = 0), the 
Whittaker function will transform into a modified, spherical Bessel function: 

W^o,ii+i/2(2«;o^) = \/ '^K'Of I T^Ki^+i/2{K,or) . In this case, the integral (8) can be 
solved exactly. 

However, we are searching for a completely analytical model which will also 
enable us to incorporate many-body effects. Our model function (9) has to be 
modified accordingly. First, we note that the asymptotic form of the Whittaker 
function as r — > 00 is 

W.,„i^^^/2i2Kor) ~ e-'^^V-"" [1 + 0{l/r)] . (10) 

Secondly, for two-body systems in which the clusters have an internal struc- 
ture, the centrifugal barrier is effectively larger and the WF should behave as 
r", where n > + 1, as r — >• 0. 

Motivated by this, we put forward the following model function 

1 6^"°'' / _ \P 

(l^b,-y{r) = -j=—^ V-~^ '^''^j ' ^^^^ 

^7=1: f^^V-ir(2«:o + m«:0'^°-'r(l-2,7^), 

m=o y^J 



(11) 



where 7 denotes the parameters {kq, Ki, r^Q,^}. The parameter kq is defined by 
the binding energy and effective mass. By putting ?7q = 770 we would ensure 
to reproduce the tail of the WF at very large r. However, higher order terms 
in Eq. (10) remain important for r < 100 fm. Therefore, r/g and are used 
as free parameters in a fit to the "exact" WF (9) in the interval of interest. 
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In this way rj^ will be an "effective" Sommerfeld parameter while ki will still 
mainly be connected with the size. Finally, the integer p is fixed by the small r 
behaviour, r", and will thus allow us to take many-body effects into account. 

It is now possible to make the integral (8) analytically (see Ref. [15]) 



V7 m=0 



X2F1 (^lf + 2 + X-r]'o;lf + l- ir]{k)] 2lf + 2; 



2ik 



rriKi + Kq + ik J 
(12) 



Many-body nuclear structure can further be taken into account by considering 
the possibility that the bound state WF contains several different two-body 
components 

Mr) = J2(^i^b,^^(^)- (13) 
i 

Note that pure many-body components will not contribute to two-body break- 
up, but will instead lead to J^if^f < 1- Note also that the threshold for two- 
body breakup will be higher for components where one (or both) of the clusters 
is excited. Therefore, we define the continuum strength function separately for 
each component. Finally, we arrive at an analytical formula for the strength 
function 



dS(EA) 



dE 



''^'(^)^^-^^^+^ ^ajk^^^^^^ (14) 



As an example we will apply our analytical approach to the ®B nucleus. Our 
key point here will be to show the applicability of our analytical model and to 
demonstrate the importance of many-body nuclear structure. The interest in 
stems from its key role in the production of high-energy solar neutrinos. The 
probability of the reaction ^Be+p — *®B+7 at solar energies strongly depends 
on the structure of ®B and, in particular, on the asymptotics of the valence 
proton WF. The reaction can be studied indirectly through Coulomb dissoci- 
ation, using a radioactive ^B beam on a heavy target [16,17,18,19]. We should 
also mention the recent progress in radiative capture measurements [20] , where 
the cross section has been measured at energies around 200 kcV with an ac- 
curacy of ~ 15 percent. However, in all cases theoretical models are needed to 
extrapolate the measured cross sections down to solar energies. 

The low-lying ^B continuum can, with relatively good precision, be approx- 
imated as a pure Coulomb one. There are no negative parity states at low 
excitation energies [21] and the electromagnetic processes are, in all cases we 
are considering, dominated by El transitions. 
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We start by treating as a pure two-body (^Be + p) system with binding 
energy Eq = 137 keV and relative orbital momentum /j = 1. The single free 
parameter, ki, in our "exact" model function (9) is then fitted to an rms 
intercluster distance of rrms = 4.57 fm (extracted from Ref. [9]). In order 
to get analytical results, we then introduce the model function (11). With 
p = 2Zj + 1 = 3 and fixed from the binding energy, the remaining parameters 
Ki and ?7o are fitted to the behavior of the "exact" model function (9), see 
Table 1. The resulting El strength function is showed as a dashed line in 
Fig. 1. This analytical approximation agrees, to a very high precision, with 
the numerical results obtained keeping the "exact" model function. The error 
is less than 2 % in the region of interest. 

However, concerning the structure of the ^B ground state, one should keep 
in mind that the ''Be core is in itself a weakly bound system with an excited 
l/2~ state at 429 keV. The common treatment of ^B as a two-body system 
is therefore highly questionable. We now want to investigate what effect the 
many-body structure of ®B might have on the strength function. We utilize 
a recent three-body calculation [9], where it was shown that, after projection 
on the two-body channel, there are three main components (adding up to 
94 % of the total WF), and that the rest are pure three-body channels, see 
Table 1. For each of these two-body components we fit our parameters ki 
and ?7q. The binding energy, Eq = 137 keV, determines kq for the two first 
components and Eq = 566 keV for the third, ^Be excited state, component. 
The best fit of the small r behavior is obtained with p = 5 which reflects the 
effectively larger centrifugal barrier in the three-body case. This centrifugal 
barrier will push the WF away from r = and will, thus, force it to become 
more narrow than the two-body WF. We therefore expect the distribution in 
momentum/energy space to be broader. This effect is clearly seen in Fig. 1 
where the El strength function obtained using this three-body modeQ is 
shown as a solid line. This difference, seen in the strength function, should 
be even more pronounced in the energy spectrum. From these results one can 
conclude that the interpretation of energy spectra is highly model dependent. 

In Fig. 2(a) we compare our analytical results for Coulomb dissociation, includ- 
ing both El and E2 transitions, to the experimental data from Davids et al. [17]. 
This experiment is very appealing since the selection of scattering angles 
{6»B < 1.77°) minimizes the contribution from nuclear scattering and the 
relatively high beam energy (82.7 MeV/A) justifies the use of first order per- 
turbation theory. Concerning the shape of the energy spectrum we have an 
excellent agreement between the experimental data and our results obtained 
using the three-body model (see thin, dotted line), while the pure two-body 



Note that this is not strictly a three-body model, but rather the two-body pro- 
jection of a three-body WF. However, in the following we will consistently refer to 
it as three-body results. 
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calculation gives a too narrow peak. As to the absolute values, the three-body 
model gives about 20 % larger cross section than the experimental data. In 
Fig. 2(b) we compare the results from our two-body model (dashed line) to 
the results from a potential model calculation without final state interaction, 
from [22] (dash-dotted). Here, the intercluster distance for both models is 
i^rms = 4.23 fm which is lower than before (rrms = 4.57 fm). The lesson from 
this figure is twofold. Firstly, it demonstrates the sensitivity to the interclus- 
ter distance and the importance of knowing this parameter. Secondly, we note 
that the difference between the two-body potential model calculation and 
our analytical results is not so large. The results deviate mainly in the peak 
amplitude. 

In connection to the same experimental conditions we show in Fig. 3 the 
fraction of the cross section that can be attributed to El transitions. This 
information is of importance when studying the inverse radiative capture re- 
action and, in contrast to the conclusion drawn by Davids et al. [17], we claim 
that although the E2 contribution increases for low relative energies, it still 
does not dominate. Note that including Ml transitions will give raise to a 
sharp dip at 640 keV. 

We have also compared our results for the radiative capture reaction to the ex- 
perimental data by Hammache et al. [20]. Our three-body result of 0-1.0(186 keV) = 
16.1 nb agrees very well with the experimental value 16.7 ±2.1 nb. We also 
find that the E2 contribution to this cross section is around 0.03 %. 

In conclusion, we have performed analytical studies of electromagnetic pro- 
cesses for loosely bound nuclei. This has been accomplished by using model 
radial functions that describe two-body WFs or the two-body projections of 
many-body WFs, accurately for all radii, and by only studying direct tran- 
sitions to/from a pure Coulomb continuum. We have examined the differ- 
ence between a pure two-body model and the results obtained incorporating 
many-body effects. From this we concluded that the interpretation of exper- 
imental data is highly model dependent. Comparisons have also been made 
to experimental results on Coulomb dissociation and, the inverse, radia- 
tive capture reaction. We found that our three-body results coincide with the 
radiative capture cross section measured by Hammache et al. [20] while our 
Coulomb dissociation energy spectrum is about 20 % larger than the experi- 
mental data from Davids et al. [17]. However, the shapes of the experimental 
and theoretical energy spectra are in excellent agreement. In contrast, our 
pure two-body model of ^B, having the same intercluster distance Trms) does 
not agree with the experimental data. A final word of wisdom is that, in 
order to interpret these data correctly, it is very important to fix the spectro- 
scopic factors of different two-body and many-body components. Therefore, 
we want to stress the usefulness of experiments where Coulomb dissociation 
is studied in complete kinematics. Examples of interesting channels in the ^B 
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case is ^B— >'''Be(l/2~)+p + 7 and ^B— >^He+a + p- Recently, '^Be fragments 
and 7-rays were measured in coincidence after breakup on a light target by 
Cortina-Gil et al. [16] and the excited core component of the WF was clearly 
observed. 

N. B. S. is grateful for support from the Royal Swedish Academy of Science. 
The support from RFBR Grants No 00-15-96590, 02-02-16174 are also ac- 
knowledged. 
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Fig. 1. The El strength function of B. Although the total strength is the same, the 
shapes are very different when treating the nucleus as a two-body (dashed line) 
or a three-body (solid line) system. The difference is mainly due to the stronger 
centrifugal barrier in the three-body case forcing the WF to be narrower and thus 
wider in momentum/energy space. The parameters of the model WFs can be found 
in Table 1. 
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Fig. 2. The ^B Coulomb dissociation energy spectrum obtained at 82.7 MeV/A 
on Pb with '^B scattering angles < 1.77°. The data points arc from [17]. In (a) 
the curves show our analytical two-body (dashed line) and three-body (solid line) 
calculations, while the thin, dotted line is just a scaled version of the three-body 
result. In (b) the dash-dotted line is the cross section calculated numerically within 
a two-body potential model, from [22]. A comparison is made with our analytical 
two-body model (dashed line) having the same intercluster distance. All theoretical 
curves have been corrected for experimental resolution and acceptance. 
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Fig. 3. Fraction of the calculated cross section for the Coulomb dissociation of 
82.7 MeV/A on Pb with scattering angles < 1.77° attributed to El transitions. 
The curves are two-body (dashed line) and three-body (solid line) calculations. 



Table 1 

Parameters of our model WFs used to describe the ground state. Both models 
give the same intercluster distance, Trms = 4.57 fm. With Eq = 137 keV we get kq = 
0.076 fm~^ and r/o = 1-595. The excited core component (last row) has Eq = 566 keV 
giving kq = 0.154 fm^^ and r/o = 0.786. The relative orbital momentum for all 
components are k = I while / is the channel spin and is the spectroscopic factor. 
Note that there is no dependence on channel spin in our strength function (14). 



Model WF 


configuration / p ni (fm ^) v'o/Vo 


two-body 


[^Be(3/2-)®Pj 


2 1.00 3 0.601 0.79 


three-body 


[7Be(3/2-)®p] 2 0.65 5 0.702 0.87 
[^Be(3/2-)®p] 1 0.13 5 0.765 0.86 
pBe(l/2-)®p] 1 0.16 5 0.753 1.43 
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